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for  the  Traveling  Salesperson  Problem  for  Dubins’  vehicle 
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Abstract — This  article  proposes  the  first  known  algorithm 
that  achieves  a  constant-factor  approximation  of  the  minimum 
length  tour  for  a  Dubins’  vehicle  through  n  points  on  the  plane. 
By  Dubins’  vehicle,  we  mean  a  vehicle  constrained  to  move  at 
constant  speed  along  paths  with  bounded  curvature  without 
reversing  direction.  For  this  version  of  the  classic  Traveling 
Salesperson  Problem,  our  algorithm  closes  the  gap  between 
previously  established  lower  and  upper  bounds;  the  achievable 
performance  is  of  order  n2/  3. 

I.  Introduction 

The  Traveling  Salesperson  Problem  (TSP)  with  its  vari¬ 
ations  is  one  of  the  most  widely  known  combinatorial 
optimization  problems.  While  extensively  studied  in  the 
literature,  these  problems  continue  to  attract  great  interest 
from  a  wide  range  of  fields,  including  Operations  Research, 
Mathematics  and  Computer  Science.  The  Euclidean  TSP 
(ETSP)  [1],  [2]  is  formulated  as  follows:  given  a  finite  point 
set  P  in  M2,  find  the  minimum-length  closed  path  through 
all  points  in  P.  It  is  quite  natural  to  formulate  this  problem 
in  context  of  Dubins’  vehicle,  i.e.,  a  non-holonomic  vehicle 
that  is  constrained  to  move  along  paths  of  bounded  curvature, 
without  reversing  direction.  The  focus  of  this  article  is  the 
analysis  of  the  TSP  for  Dubins’  vehicle;  we  shall  refer  to  it 
as  DTSP. 

Exact  algorithms,  heuristics  as  well  as  polynomial-time 
constant  factor  approximation  algorithms  are  available  for 
the  Euclidean  TSP,  see  [3],  [4],  [5].  It  is  known  that  non¬ 
metric  versions  of  the  TSP  are,  in  general,  not  approximable 
in  polynomial  time  [6],  Furthermore,  unlike  most  other 
variations  of  the  TSP,  it  is  believed  that  the  DTSP  cannot  be 
formulated  as  a  problem  on  a  finite-dimensional  graph,  thus 
preventing  the  use  of  well-established  tools  in  combinatorial 
optimization.  On  the  other  hand,  it  is  reasonable  to  expect 
that  exploiting  the  geometric  structure  of  Dubins’  paths  one 
can  gain  insight  into  the  nature  of  the  solution,  and  possibly 
provide  polynomial-time  approximation  algorithms. 

A  fairly  complete  picture  is  available  for  the  minimum¬ 
time  point-to-point  path  planning  problem  for  Dubins’  ve¬ 
hicle,  see  [7]  and  [8].  However,  the  DTSP  seems  not  to 
have  been  studied  as  extensively.  In  [9],  some  results  for  the 
worst  case  tours  of  DTSP  were  provided.  A  lower  bound 
on  the  expected  cost  of  a  stochastic  DTSP  visiting  ran¬ 
domly  generated  points  was  provided  in  [10].  Here,  we  shall 
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specifically  concentrate  on  the  case  when  the  target  points 
in  the  environment  are  generated  stochastically  according 
to  a  uniform  distribution.  We  shall  refer  to  such  a  problem 
as  stochastic  DTSP.  In  this  context  the  first  algorithm  with 
asymptotic  sub-linear  cost  was  proposed  in  [1 1];  an  improved 
algorithm  was  proposed  in  [12]. 

The  motivation  to  study  the  DTSP  arises  in  robotics 
and  uninhabited  aerial  vehicles  (UAVs)  applications,  e.g., 
see  [13],  [14],  [15].  In  particular,  we  envision  applying  our 
algorithm  to  the  setting  of  an  UAV  monitoring  a  collection  of 
spatially  distributed  points  of  interest.  Additionally,  from  a 
purely  scientific  viewpoint,  it  appears  to  be  of  general  interest 
to  bring  together  the  work  on  Dubins’  vehicle  and  that  on 
TSP.  Some  concrete  results  along  these  lines  have  been 
obtained  in  [9]  and  in  [11]  where  an  algorithm  to  guarantee 
sub-linear  cost  for  the  stochastic  DTSP  was  proposed.  UAV 
applications  also  motivate  us  to  study  the  Dynamic  Traveling 
Repairperson  Problem  (DTRP),  in  which  the  aerial  vehicle 
is  required  to  visit  a  dynamically  changing  set  of  targets. 
This  problem  was  introduced  by  Bertsimas  and  van  Ryzin 
in  [16]  and  then  decentralized  policies  achieving  the  same 
performances  were  proposed  in  [13].  However,  as  with  the 
TSP,  the  study  of  DTRP  in  context  of  Dubins’  vehicle  has 
eluded  attention  from  the  research  community. 

The  contributions  of  this  article  are  twofold.  First,  we 
propose  an  algorithm  for  the  stochastic  DTSP  through  a  point 
set  P,  called  the  Recursive  Bead-Tiling  Algorithm, 
based  on  a  geometric  tiling  of  the  plane,  tuned  to  the 
Dubins’  vehicle  dynamics,  and  a  strategy  for  the  vehicle 
to  service  targets  from  each  tile.  Second,  we  obtain  an 
upper  bound  on  the  stochastic  performance  of  the  proposed 
algorithm  and  thus  also  establish  a  similar  bound  on  the 
stochastic  DTSP.  The  upper  bound  on  the  performance 
of  the  Recursive  Bead-Tiling  Algorithm  belongs  to 
0(n2/3)  with  high  probability,  and  it  is  known  that  the 
lower  bound  on  the  achievable  performance  belongs  to 
Q  (n2/3).  The  algorithm  we  introduce  in  this  article  is  the 
first  known  algorithm  providing  a  provable  constant-factor 
approximation  to  the  DTSP  optimal  solution. 

Notation 

Here  we  collect  some  concepts  that  will  be  required  in 
the  later  sections.  A  Dubins’  vehicle  is  a  planar  vehicle  that 
is  constrained  to  move  along  paths  of  bounded  curvature, 
without  reversing  direction  and  maintaining  a  constant  speed. 
Accordingly,  we  define  a  feasible  curve  for  Dubins’  vehicle 
or  a  Dubins’  path,  as  a  curve  that  is  twice  differentiable 
almost  everywhere,  and  such  that  the  magnitude  of  its 


curvature  is  bounded  above  by  1/p,  where  p  >  0  is  the 
minimum  turn  radius. 

Let  P  =  {pi, . . .  ,pn}  be  a  set  of  n  points  in  a  compact 
region  Q  C  R2  and  P„  be  the  collection  of  all  point  sets 
P  C  Q  with  cardinality  n.  Let  ETSP(P)  denote  the  cost  of 
the  Euclidean  TSP  over  P,  i.e.,  the  length  of  the  shortest 
closed  path  through  all  points  in  P.  Correspondingly,  let 
DTSPp(P)  denote  the  cost  of  the  Dubins’  TSP  over  P,  i.e., 
the  length  of  the  shortest  closed  Dubins’  path  through  all 
points  in  P.  In  what  follows,  p  £  R+  is  take  constant,  and 
we  study  the  dependence  of  DTSPp  :  Vn  — >  R+  on  n. 

For  /,  g  :  N  — >  R,  we  say  that  /  £  0(g)  (respectively, 
/  £  f2(</))  if  there  exist  Nq  £  N  and  k  £  R+  such  that 
| /(iV) |  <  k\g(N)\  for  all  N  >  N0  (respectively,  \f(N)\  > 
k\g(N)\  for  all  N  >  Nq).  If  f  £  0(g)  and  /  £  f 1(g),  then 
we  use  the  notation  /  £  0(g). 

II.  The  stochastic  DTSP 

In  [9],  a  simple  heuristics,  the  ALTERNATING  ALGO¬ 
RITHM  for  the  Dubins’  TSP  for  a  given  point  set  was 
proposed.  The  length  of  tour  generated  by  this  algorithm 
was  also  characterized  and  it  was  shown  that  it  belongs  to 
fl(^/n)  and  0(n).  It  was  also  shown  that  this  simple  policy 
performs  well  when  the  points  to  be  visited  by  the  tour  are 
chosen  in  an  adversarial  manner.  However,  it  is  reasonable  to 
argue  that  this  algorithm  might  not  perform  very  well  when 
dealing  with  a  random  distribution  of  the  target  points.  In 
particular,  one  can  expect  that  when  n  points  are  chosen 
randomly,  the  cost  of  the  DTSP  increases  sub-linearly  with 
n,  i.e.,  that  the  average  length  of  the  path  between  two  points 
decreases  as  n  increases.  In  this  section,  we  consider  the 
scenario  when  n  target  points  are  stochastically  generated  in 
Q  according  to  a  uniform  distribution.  A  novel  algorithm, 
the  Bead-Tiling  Algorithm  was  proposed  in  [11]  to 
service  these  points  in  such  a  way  that  its  tour  length  grew 
sub-linearly  with  the  number  of  points  asymptotically  with 
high  probability,  where  an  event  is  said  to  occur  with  high 
probability  if  the  probability  of  its  occurence  approaches  1  as 
n  — >  Too.  Here,  we  present  a  novel  version  of  this  strategy 
in  the  form  of  the  Recursive  Bead-Tiling  Algorithm 
and  characterize  its  performance. 

We  assume  that  the  environment  Q  is  a  rectangle  of  width 
W  and  height  // ;  different  choices  for  the  shape  of  Q  affect 
our  conclusions  only  by  a  constant.  In  what  follows  we  select 
a  reference  frame  whose  two  axes  are  parallel  to  the  sides  of 
Q.  Let  n  target  points  be  generated  stochastically  according 
to  uniform  distribution  in  the  region  Q.  Let  V  =  (pi, . . .  ,pn) 
be  the  locations  of  these  target  points. 

A.  A  lower  bound 

First,  we  summarize  a  result  from  [10],  that  provides  a 
lower  bound  on  the  expected  length  of  the  stochastic  DTSP. 

Theorem  2.1:  (Lower  bound  on  stochastic  DTSP)  For  all 
p  >  0,  the  expected  cost  of  a  stochastic  DTSP  visiting 
a  set  P  of  n  uniformly-randomly-generated  points  in  Q, 
E[DTSPp(P)]  belongs  to  f 2(n2/3). 


B.  A  constructive  upper  bound 

In  [11],  a  novel  algorithm,  the  Bead-Tiling  Algo¬ 
rithm,  to  compute  Dubins’  path  through  a  point  set  in 
the  region  Q  was  proposed.  In  this  section,  we  design  the 
Recursive  Bead-Tiling  Algorithm  based  on  the  ideas 
of  the  Bead-Tiling  Algorithm.  We  will  show  that  the 
proposed  algorithm  provides  a  constant  factor  approximation 
to  the  optimal  DTSP  with  high  probability.  We  start  by 
describing  some  useful  geometric  objects. 

1 )  The  basic  geometric  construction:  Consider  two  points 
=  (— /,  0)  and  p+  =  (l,  0)  on  the  plane,  with  l  <  p, 
and  construct  the  region  Bp(l)  as  detailed  in  Figure  |T|  In 


Fig.  1.  Construction  of  the  “bead”  The  figure  shows  how  the  upper 

half  of  the  boundary  is  constructed,  the  bottom  half  is  symmetric. 


the  following,  we  will  refer  to  such  regions  as  beads.  The 
region  Bp(l)  enjoys  the  following  asymptotic  properties  as 
(l/p)  -  0+: 

(PI)  The  maximum  “thickness”  of  the  region  is  equal  to 

(P2)  The  area  of  Bp(l)  is  equal  to 

Area[Sp(Z)]  =  lw(l)  =  l—  T  o  j  . 

(P3)  For  any  p  £  Bp,  there  is  at  least  one  Dubins’  path 
7 p  through  the  points  {p-,p,p+},  entirely  contained 
within  Bp,  and  such  that  its  length  is  at  most 

Length(7p)  <  4parcsin  ^  =  21  +  o 

These  facts  are  verified  using  elementary  planar  geometry. 

2 )  Periodic  tiling  of  the  plane:  An  additional  property  of 
the  geometric  shape  introduced  above  is  that  the  plane  can 
be  periodically  tiled  by  identical  copies  of  Bp(l),  for  any 
l  £  (0,2 p\.  (Recall  that  a  tiling  of  the  plane  is  a  collection 


of  sets  whose  intersection  has  measure  zero  and  whose  union 
covers  the  plane.)  Let 


_  Arca[gp(Q] 

M  j  Area[Q]  ' 

Consider  a  bead  B  entirely  contained  in  Q;  the  probability 
that  the  f-th  point  is  sampled  in  B  is  equal  to  //.  Furthermore, 
the  probability  that  exactly  k  out  of  the  n  points  are  sampled 
in  B  has  a  binomial  distribution,  i.e.,  indicating  with  ns  the 
total  number  of  points  sampled  in  B, 


Prob[riB  =  k\n  samples] 


(1) 


Choose  /x  as  a  function  of  n,  in  such  a  way  that  v  =  n/j,(n) 
is  a  constant.  In  such  a  case,  the  limit  for  large  n  of  the 
binomial  distribution  Q  is  the  Poisson  distribution  of  mean 
v,  that  is, 

l jk 

lim  Prob[ns  =  k\n  samples]  =  — e  "  (2) 

n— >oo  h\ 

C.  The  algorithm 

Consider  a  tiling  of  the  plane  such  that  Area[£>p(7)]  = 
WH/{2n)\  in  such  a  case,  /i  =  l/(2n),  and  v  =  1/2.  (Note 
that  this  implies  that  n  must  be  large  enough  that  l  <  2 p.) 
Furthermore,  the  tiling  is  chosen  is  such  a  way  that  it  is 
aligned  with  the  sides  of  Q,  see  Figure  [2] 


Fig.  2.  Sketch  of  tiling  of  the  region  before  the  first  phase  of  the 
Recursive  Bead-Tiling  Algorithm. 

The  proposed  algorithm  will  consist  of  a  sequence  of 
phases;  during  each  of  these  phases,  a  Dubins  tour  (i.e.,  a 
closed  path  with  bounded  curvature)  will  be  constructed  that 
“sweeps”  the  set  Q. 

In  the  first  phase,  a  Dubins  tour  is  constructed  with  the 
following  properties: 

(i)  it  visits  all  non-empty  beads  once, 

(ii)  it  visits  all  rows1  in  sequence  top-to-down,  alternat¬ 
ing  between  left-to-right  and  right-to-left  passes,  and 
visiting  all  non-empty  beads  in  a  row, 

(iii)  when  visiting  a  non-empty  bead,  it  services  at  least 
one  target  in  it. 

In  order  to  visit  the  outstanding  targets,  a  new  phase  is 
initiated.  In  this  phase,  instead  of  considering  single  beads, 

*A  row  is  a  maximal  string  of  beads  with  non-empty  intersection  with 

<2. 


we  will  consider  “meta-beads”  composed  of  two  beads  each, 
as  shown  in  Figure  m,  and  proceed  in  a  similar  way  as 
the  first  phase,  i.e.,  a  Dubins  tour  is  constructed  with  the 
following  properties: 

(i)  the  tour  visits  all  non-empty  meta-beads  once, 

(ii)  it  visits  all  (meta-bead)  rows  in  sequence  top-to- 
down,  alternating  between  left-to-right  and  right-to-left 
passes,  and  visiting  all  non-empty  meta-beads  in  a  row, 

(iii)  when  visiting  a  non-empty  meta-bead,  it  services  at 
least  one  target  in  it. 

This  process  is  iterated  at  most  log2  n  +  1  times,  and 
at  each  phase  meta-beads  composed  of  two  neighboring 
meta-beads  from  the  previous  phase  are  considered;  in  other 
words,  the  meta-beads  at  the  f-th  phase  are  composed  of  2l~1 
neighboring  beads.  After  the  last  phase,  the  leftover  targets 
will  be  visited  using,  for  example,  a  greedy  strategy,  or  the 
Alternating  Algorithm. 

We  have  the  following  result,  which  we  prove  using  a 
technique  similar  to  that  developed  in  [17]. 

Lemma  2.2:  Let  P  £  Vn  be  uniformly  randomly  gener¬ 
ated  in  Q.  Then,  the  number  of  unvisited  targets  after  the 
last  phase  of  the  Recursive  Bead-Tiling  Algorithm 
over  P  belongs  to  O(logro)  with  high  probability. 

Proof:  Associate  a  unique  identifier  to  each  bead,  e.g., 
integers  between  1  and  2 n;  call  such  a  set  of  identifiers  X.  Let 
b(t )  €  X  be  the  identifier  of  the  bead  in  which  the  f-th  target 
is  sampled,  and  let  h(t)  £  N  be  the  phase  at  which  the  f-th 
target  is  visited.  Without  loss  of  generality,  we  will  assume 
that  if  b(t\)  =  biff),  and  t\  <  f2,  then  h{t{)  <  h{tf). 

Indicate  with  v-ff)  the  number  of  beads  that  contain  unvis¬ 
ited  targets  at  the  inception  of  the  f- th  phase,  computed  after 
the  insertion  of  the  f-th  target.  Furthermore,  let  rrii  be  the 
number  of  f-th  phase  meta-beads  (i.e.,  meta-beads  containing 
2l_1  neighboring  beads)  with  a  non-empty  intersection  with 
Q.  Clearly,  Vi(t)  <  Uj(n),  rrii  <  2mi+i,  and  v\ (n)  <  n  < 
mi/ 2  with  certainty. 

The  f-th  target  will  not  be  visited  during  the  first  phase  if 
it  is  sampled  in  a  bead  that  already  contains  other  targets.  In 
other  words. 


Pr[fi(f)  >  2|ui(f)] 


Mt)  K  vi(p)  <  l 

mi  ~  2n  ~  2 


Similarly,  the  f-th  target  will  not  be  visited  during  the  f-th 
phase  if  (i)  it  has  not  been  visited  before  the  f-th  pass,  and  (ii) 
it  belongs  to  a  meta-bead  that  already  contains  other  targets 
not  visited  before  the  f-th  phase: 


Pr  [h(t)  >  i  +  l\{vi(t  -  l),^_i(f  -  l),Ui(f  -  1))] 
=  Pr[/i(f)  >  f  +  l\h(t)  >  f,  Vi(t  —  1)]  • 
Pr[/i(f)  >  i\(vi-i(t-  1),  ■  •  • ,  Vi(t  —  1))] 
Vi(t  -  1) 


< 


n 

3  =  1 


rrii 

rrij 


■Pr[ft(f)  >  f |(v*_i(f  —  1), . . .  ,vi(t  -  1))] 


< 


-q  2-?  1vj(n) 


'2V 


3= 1 


2n 


n  <*(»)■ 

3= 1 


(3) 


Phase  1 


Phase  2 


Phase  3 


Fig.  3.  Sketch  of  “meta-beads”  at  successive  phases  in  the  recursive  bead  tiling  algorithm. 


For  a  fixed  i  >  1,  define  a  sequence  of  binary  random 
variables 

Y  _  /  1)  if  h(t)  >  i  +  1  and  vfat  -  1)  <  fan] 

4  (  0,  otherwise. 

In  other  words,  Yt  is  equal  to  1  if  the  /-th  target  is  not  visited 
within  the  first  i  phases  despite  the  fact  that  the  number  of 
beads  still  containing  unvisited  target  at  the  inception  of  the 
7-th  phase  is  less  than  /3jn;  the  values  {A}  will  be  defined 
shortly. 

Even  though  the  random  variable  Yt  depends  on  the  targets 
generated  before  the  t-th  target,  the  probability  that  it  takes 
the  value  1  is  bounded  by 


Pr [Yt  =  1|6(1), 6(2),  •  •  • ,  b(t  -  1)]  <  2^  H  fa  =:  Pi, 

3= 1 


regardless  of  the  actual  values  of  6(1), . . .  ,b(t  —  1). 

It  is  known  (e.g.,  see  [17])  that  if  the  random  variables 
Yt  satisfy  such  a  condition,  the  sum  Yt  is  stochastically 
dominated  by  a  binomially  distributed  random  variable, 
namely, 


Pr 


J2Yt>k 

_t= l 


<  Pr [B(n,pi)  >  fc]. 


In  particular. 


which  is  less  than  1/n2  for  i  <  i*(n )  :=  [l°g2 n  — 
log2  log2  n  —  log2  6J  <  log2  n.  Note  that 

A  <12^  V*  >  i*(n).  (7) 

n 

Let  £j  be  defined  as  the  event  that  Vi(n)  <  fan.  Note  that 
if  A  is  true,  then  Ui+i(n)  <  Yt:  the  right  hand  side 

represents  the  number  of  targets  that  will  be  visited  after  the 
7-th  phase,  whereas  the  left  hand  side  counts  the  number  of 
beads  containing  such  targets.  We  have,  for  all  i  <  i*(n ): 


Pr[i>i+i  >  /3i+H£i]’Pr[£i]  <  Pr 


^Yt>  fa+i n 

_t=i 


1 

<  , 


that  is  , 


Pr  {-^Si+i\Si}  < 


1 

n2  Pr[fi]  ’ 


and  thus  (recall  that  S\  is  true  with  certainty): 


Pr  [-i£i+i]  <  \  +  Pr[-i£i]  < 

n  nz 

In  other  words,  for  all  i  <  i*(n),  vfan)  <  fan  with  high 
probability. 

Let  us  turn  our  attention  to  the  phases  such  that  i  >  i*(n). 
The  total  number  of  targets  visited  after  the  i*- th  phase 
is  dominated  by  a  binomial  variable  B(n,  12  log2  n/n)\  in 
particular, 


Pr 


^  Yt  >  2nPi 


<  Pr [B(n,pi)  >  2npi \  <  2  npi/3, 


(5) 

where  the  last  inequality  is  obtained  using  Chernoff’s  bound. 
Let  us  define  the  sequence  {fa}  through  the  recursion 


fa  —  1) 


Hi- 3) 


fa+i  =  2Pl  =  2^+1\{fa  =  ?-2  pi 

3=1 

which  leads  to 


A  =  2 


1-2 


(6) 


With  the  above  definition  in  mind,  (0  can  be  rewritten 


as 


Pr 


y:  Yt  >  Pi+ii 


,4=1 


<  Pr [B(n,pi]  >  fa+in] 

C  2~Pi+in/6  =  2~T¥ 


Pr  [ttj.+i  >  24  log2  n|A»]-Pr[A*]  <  Pr 


^  Yf  >  24  log 2n 
.4=1 


<  Pi'[B(n,  12  log2  n/n)  >  241og2n]  <  2  121og2ra; 


dealing  with  conditioning  as  before,  we  get 

PrIW  >  24 ^ 2"]  <  ^J+PrKi*]  <  ^  +  ^|^- 

In  other  words,  the  number  of  targets  that  will  be  left  after 
the  **-th  phase  will  be  bounded  by  a  logarithmic  function  of 
n  with  high  probability.  ■ 


D.  A  bound  on  the  length  of  the  solution 

What  we  know  at  this  point  is  that  after  a  sufficiently 
large  number  of  phases,  almost  all  targets  will  be  visited, 
with  high  probability.  The  key  point  is  to  recognize  that  the 
length  of  each  phase  is  decreasing  at  such  a  rate  that  the  sum 
of  the  lengths  of  all  the  phases  remains  bounded.  We  first 
state  and  prove  the  following  result  which  characterizes  the 


length  of  Dubins’  path  required  to  execute  the  RECURSIVE 
Bead-Tiling  Algorithm. 

Lemma  2.3:  (Length  of  path  for  the  Recursive  Bead- 
Tiling  Algorithm)  Let  P  e  Vn  be  uniformly  randomly 
generated  in  Q.  Then  the  length  of  Dubins’  path  required 
to  execute  logn  phases  of  the  RECURSIVE  Bead-Tiling 
Algorithm  over  P  belongs  to  0(n2/3). 

Proof:  Let  L,  denote  the  upper  bound  on  the  length  of 
the  path  for  the  i*  phase.  Then  one  can  see  that 

log2(») 

L,<  3  E  L2j-i- 
i= 1  j=l 

Let  us  first  compute  the  length  of  a  pass,  in  either 
direction.  The  number  of  beads  traversed  will  be  no  more 
than 


r  w  i 

Ci  il 

2Hn 

— n3 
2i 

where  c\  =  i/fyn  is  a  constant.  The  length  of  Dubins’ 
path  contained  entirely  within  a  meta-bead  at  the  (2 j  —  l)-th 
phase  is  less  than  2-7-1  (2 ln  +  o((2)).  Hence,  the  total  path 
length  per  pass  will  be  bounded  by: 

Tpass,2j— l  <  1  (2(n  +  o((2))  ro3  +  1^) 

=  c\lnn^  +  2Jln  +  ^-ri3o((2)  +  2-J_1o((2) 

(9) 


as  !„  ->  0+. 

The  cost  of  a  U-turn,  i.e.,  the  length  of  the  path  needed 
to  reverse  direction  and  move  to  the  next  row  of  beads,  is 
bounded  by 


Pu-turn,2j-l  <  ^p+2J  2  w{ln)  =  ^-TT  p+23  2(^f-+o((3)). 

(10) 

The  total  number  of  passes  will  be  at  most 


Nx 


pass,2j  — 1 


H 


2-1  2w(ln) 


< 


pH 

2 


+  1. 


(11) 


The  cost  of  closing  the  tour  is  bounded  by  a  constant,  say 


-^closure, 2j—  1  <  4(W  +  Http).  (12) 


Concluding,  the  total  path  length  will  be  bounded  by 


L‘2j  —  1  —  ^pass,2j  — 1  (Ppass^j  —  lPPu— turn, 2j  —  l)“l_ -^closure, 2j-l- 

(13) 

Substituting  eqns.  <!()>.  (lilt  and  (1121)  in  eqn.  j  1 3t.  one 
can  find  constants  k\ ,  F'2  and  k3  such  that 


L2j-i  <  ki(l+n^o(n~i)+n^)+k22j  +o(n"5)j 

+  k32~j  (l  +  no(n~ 3 )  +  n3) . 

From  this  expression  for  the  length  of  path  during  odd 
phases,  one  can  conclude  that 

iog(n)  r  log2(,1)  1 

E  P;  <  3  E  L2j- i  G  0(n2/3). 

i=i  j= l 


Based  on  the  results  obtained  so  far,  we  are  now  ready  to 
state  an  upper  bound  on  the  length  of  the  path  traveled  by 
Dubins’  vehicle  to  service  all  the  n  targets  while  executing 
the  Recursive  Bead-Tiling  Algorithm  followed  by 
the  Alternating  Algorithm;  let  LReCBTA,P(P)  repre¬ 
sent  the  corresponding  quantity. 

Theorem  2.4:  (Upper  bound  on  the  length  of  the  total 
path)  Let  F  £  P„  be  uniformly  randomly  generated  in 
Q.  For  all  p  >  0,  LRecBTA,P(P)  G  0(n2/3)  with  high 
probability. 

Proof:  By  Lemma  12.31  the  length  of  path  to  exe¬ 
cute  the  Recursive  Bead-Tiling  Algorithm  belongs 
to  0(n2/3).  From  Theorem  12.21  the  number  of  targets 
remaining  at  the  end  of  the  Recursive  Bead-Tiling 
Algorithm  belongs  to  O(logn)  with  high  probability. 
These  remaining  points  can  be  serviced  by  any  greedy  policy 
or  some  heuristics  (e.g.,  Alternating  Algorithm  [9]) 
in  O(logn)  time.  The  statement  of  the  theorem  follows 
immediately.  ■ 

Combining  results  from  Theorem  12.11  and  Theorem  12.41 
one  can  conclude  that  the  Recursive  Bead-Tiling  Al¬ 
gorithm  is  a  constant  factor  approximation  to  the  optimal 
DTSP  with  high  probability. 

III.  Conclusions 

In  this  article,  we  have  studied  the  TSP  problem  for  vehi¬ 
cles  that  follow  paths  of  bounded  curvature  in  the  plane.  For 
the  stochastic  setting,  we  have  obtained  upper  bounds  that 
are  within  a  constant  factor  of  the  lower  bound  established 
in  literature  [10];  the  upper  bounds  are  constructive  in  the 
sense  that  they  are  achieved  by  two  novel  algorithms.  It 
is  interesting  to  compare  our  results  with  the  Euclidean 
setting  (i.e.,  the  setting  in  which  curves  do  not  have  curvature 
constraints).  For  a  given  compact  set  and  a  point  set  P  of 
n  points,  it  is  known  [1],  [2]  that  the  ETSP(P)  belongs 
to  Q(y/n).  This  is  true  for  both  stochastic  and  worst-case 
settings.  In  this  article,  we  showed  that,  given  a  fixed  p  >  0, 
the  stochastic  DTSPp(P)  belongs  to  @(n2/3)  with  high 
probability.  It  is  known  [9]  that  the  worst-case  DTSPp(P) 
belongs  to  0(n). 

In  the  future,  we  plan  to  perform  extensive  simulations  to 
support  the  results  obtained  in  this  article.  Future  directions 
of  research  include  study  of  centralized  and  decentralized 
versions  of  the  DTRP  and  general  task  assignment  and 
surveillance  problems  for  various  non-holonomic  vehicles. 

Acknowledgment 

This  material  is  based  upon  work  supported  in  part  by 
ONR  YIP  Award  N000 14-03 -1-05 12  and  AFOSR  MURI 
Award  F49620-02- 1-0325.  The  authors  would  like  to  thank 
John  J.  Enright  for  helpful  discussions. 

References 

[1]  J.  Beardwood,  J.  Halton,  and  J.  Hammersly,  “The  shortest  path  through 
many  points,”  in  Proceedings  of  the  Cambridge  Philosophy  Society, 
vol.  55,  pp.  299-327,  1959. 


[2]  J.  M.  Steele,  “Probabilistic  and  worst  case  analyses  of  classical  prob¬ 
lems  of  combinatorial  optimization  in  Euclidean  space,”  Mathematics 
of  Operations  Research ,  vol.  15,  no.  4,  p.  749,  1990. 

[3]  D.  Applegate,  R.  Bixby,  V.  Chvatal,  and  W.  Cook,  “On  the  solu¬ 
tion  of  traveling  salesman  problems,”  in  Documenta  Mathematica, 
Journal  der  Deutschen  Mathematiker-Vereinigung,  (Berlin,  Germany), 
pp.  645-656,  Aug.  1998.  Proceedings  of  the  International  Congress 
of  Mathematicians,  Extra  Volume  ICM  III. 

[4]  S.  Arora,  “Nearly  linear  time  approximation  scheme  for  Euclidean  TSP 
and  other  geometric  problems,”  in  Proc.  38th  IEEE  Annual  Symposium 
on  Foundations  of  Computer  Science,  (Miami  Beach,  FL),  pp.  554- 
563,  Oct.  1997. 

[5]  S.  Lin  and  B.  W.  Kemighan,  “An  effective  heuristic  algorithm  for  the 
traveling-salesman  problem,”  Operations  Research,  vol.  21,  pp.  498- 
516,  1973. 

[6]  S.  Sahni  and  T.  Gonzalez,  “P-complete  approximation  problems,” 
Journal  of  the  Association  of  Computing  Machinery,  vol.  23,  no.  3, 
pp.  555-565,  1976. 

[7]  L.  E.  Dubins,  “On  curves  of  minimal  length  with  a  constraint  on 
average  curvature  and  with  prescribed  initial  and  terminal  positions 
and  tangents,”  American  Journal  of  Mathematics,  vol.  79,  pp.  497- 
516,  1957. 

[8]  A.  M.  Shkel  and  V.  J.  Lumelsky,  “Classification  of  the  Dubins  set,” 
Robotics  and  Autonomous  Systems,  vol.  34,  pp.  179-202,  2001. 

[9]  K.  Savla,  E.  Frazzoli,  and  F.  Bullo,  “On  the  point-to-point  and 
traveling  salesperson  problems  for  Dubins’  vehicle,”  in  American 
Control  Conference,  (Portland,  OR),  pp.  786-791,  June  2005. 

[10]  J.  J.  Enright  and  E.  Frazzoli,  “UAV  routing  in  a  stochastic  time-varying 
environment,”  in  IFAC  World  Congress,  (Prague,  Czech  Republic), 
July  2005.  Electronic  Proceedings. 

[11]  K.  Savla,  F.  Bullo,  and  E.  Frazzoli,  “On  traveling  salesperson  problems 
for  Dubins’  vehicle:  stochastic  and  dynamic  environments,”  in  IEEE 
Conf  on  Decision  and  Control,  (Seville,  Spain),  pp.  4530^1-535,  Dec. 
2005. 

[12]  S.  Itani  and  M.  A.  Dahleh,  “On  the  average  length  of  the  stochastic 
TSP  with  Dubin’s  metric.”  Submitted  to  the  American  Control  Con¬ 
ference,  2005. 

[13]  E.  Frazzoli  and  F.  Bullo,  “Decentralized  algorithms  for  vehicle  routing 
in  a  stochastic  time-varying  environment,”  in  IEEE  Conf  on  Decision 
and  Control,  (Paradise  Island,  Bahamas),  pp.  3357-3363,  Dec.  2004. 

[14]  R.  W.  Beard,  T.  W.  McLain,  M.  A.  Goodrich,  and  E.  P.  Anderson, 
“Coordinated  target  assignment  and  intercept  for  unmanned  air  vehi¬ 
cles,”  IEEE  Transactions  on  Robotics  and  Automation,  vol.  18,  no.  9, 
pp.  911-922,  2002. 

[15]  S.  Darbha,  “Combinatorial  motion  planning  for  a  collection  of  Reeds- 
Shepp  vehicles,”  tech,  rep.,  ASEE/AFOSR  SFFP,  AFRL,  Eglin,  Aug. 
2005. 

[16]  D.  J.  Bertsimas  and  G.  J.  van  Ryzin,  “A  stochastic  and  dynamic  vehicle 
routing  problem  in  the  Euclidean  plane,”  Operations  Research,  vol.  39, 
pp.  601-615,  1991. 

[17]  Y.  Azar,  A.  Z.  Broder,  A.  R.  Karlin,  and  E.  Upfal,  “Balanced 
allocations,”  SIAM  Journal  on  Computing,  vol.  29,  no.  1,  pp.  1 80— 
200,  1999. 


